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New predictions for the anti-neutrino flux emitted by nuclear reactors suggest that reactor ex- 
periments may have measured a deficit in the anti-neutrino flux, which can be interpreted in terms 
of oscillations between the known active neutrinos and new sterile states. Motivated by this obser- 
vation, we perform a re-analysis of global short-baseline neutrino oscillation data in a framework 
with one or two sterile neutrinos. While one sterile neutrino is still not sufficient to reconcile the 
signals suggested by reactor experiments and by the LSND and MiniBooNE experiments with null 
results from other searches, we find that, with the new reactor flux prediction, the global fit improves 
considerably when the existence of two sterile neutrinos is assumed. 

PACS numbers: 14.60.Pq, 14.60.St 
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Introduction. By now a standard paradigm of neutrino 
physics has emerged. A beautiful series of experiments 
has established the phenomenon of neutrino oscillations. 
Results from solar, atmospheric, reactor, and accelerator 
neutrino experiments can be accommodated nicely by os- 
cillations of the three neutrinos of the Standard Model, 
the so-called "active" neutrinos, with mass-squared dif- 
ferences of order 10~^ and 10~^ eV^, see [H, for recent 
fits and references. There are, however, a few exper- 
imental results which cannot be explained within this 
framework and seem to require additional neutrinos with 
masses at the eV scale 0, 01 ■ Such neutrinos cannot 
participate in the weak interactions due to collider con- 
straints, and are therefore called "sterile" neutrinos. 

Recently another hint for sterile neutrinos has emerged 
from a re-evaluation of the expected anti-neutrino flux 
emitted from nuclear reactors [5|. The new prediction is 
about 3% higher than what was previously assumed Q. 
If confirmed, this result would imply that all existing 
neutrino oscillation searches at nuclear reactors have ob- 
served a deficit of electron anti- neutrinos (Pg), which can 
be interpreted in terms of oscillations at baselines of or- 
der 10-100 m 0. At typical reactor anti-neutrino ener- 
gies of few MeV, standard oscillations of the three active 
neutrinos require baselines of a least 1 km. Hence, the 
"reactor anomaly" can only be accommodated if at least 
one sterile neutrino with mass at the eV-scale or higher is 
introduced. This is particularly intriguing because also 
the long-standing "LSND anomaly" [3|, as well as the 
more recent MiniBooNE anti-neutrino results Q suggest 
the existence of a sterile neutrino in that mass range. 

Previous phenomenological studies have been 

performed in a framework in which the standard three ac- 
tive neutrino scenario is amended by adding one ("3-1-1") 
or two ("3-1-2") sterile neutrinos with masses at the eV 
scale. These studies came to the conclusion that an ex- 
planation of the aforementioned anomalies within these 
sterile neutrino scenarios is in conflict with various con- 
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xVdof 


3+1 


1.78 


0.151 




50.1/67 


3+2 


0.46 


0.108 0.89 


0.124 


46.5/65 



Table I: Best fit points for the 3+1 and 3+2 scenarios from 
reactor anti-neutrino data. The total number of data points 
is 69 (Bugey3 spectra plus 9 SBL rate measurements; we have 
omitted data from Chooz and Palo Verde, which are not very 
sensitive to the model parameters, but would dilute the by 
introducing 15 additional data points). For no oscillations we 
have xVdof = 59.0/69. 



straints from other neutrino oscillation searches at short 
baselines (SBL), including also data from reactor exper- 
iments. In this note we revisit 3+1 and 3+2 sterile neu- 
trino oscillation schemes in the light of the new reactor 
neutrino fluxes. We argue that one sterile neutrino is still 
not sufficient to describe all data, whereas a 3+2 frame- 
work is now in much better agreement with the data. 

A^ew reactor fluxes and fit of SBL reactors. Let us 
first discuss the implications of the new reactor anti- 
neutrino flux prediction for reactor data alone by an- 
alyzing a set of SBL reactor experiments at baselines 
L < 100 m Q. We include full spectral data from the 
Bugcy3 experiment 11 1| at 15, 40 and 95 m and take into 
account the Bugey4 ]i1l, ROVNO [H, Krasnoyarsk [li], 
ILL [l5| . and Gosgen |l6| experiments via the rate mea- 
surements summarized in Table II of Q- Furthermore 
we include the Chooz 17 1 and Palo Verde [l^l experi- 
ments at L ~ 1 km. We use the neutrino fluxes from the 



isotopes 235u, 239p^^ 238u^ 24ipy obtained in ^ and we 
include the uncertainty on the integrated flux for each 
isotope given in Table I of , correlated between all ex- 
periments. For further technical details see [ij. 

We perform a flt to these data within the 3+1 and 
3+2 sterile neutrino frameworks, where neutrino oscil- 
lations for SBL reactor experiments depend on 2 and 4 
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Figure 1: Comparison of sterile neutrino models to reactor 
data: energy spectra from BugeyS and the rate measurement 
of Bugey4 (inset). The data points correspond to the ratio 
of the observed event numbers to the predicted event num- 
ber for no oscillations using the new reactor anti-neutrino 
fluxes For the BugeyS spectra we show statistical er- 

rors only, whereas the error on the Bugey4 rate is dominated 
by systematics. The green solid curve shows the prediction 
for the no oscillation hypothesis, the blue solid and red solid 
curves correspond to the 3-1-1 and 3-1-2 best fit points for SBL 
reactor data (Tab. [J), and the dashed curve corresponds to 
the 3-1-2 best fit point of global SBL data (Tab. Hill. 



parameters, respectively. The parameters are the mass- 
squared differences Am^^ and Am|]^ between the eV- 
scale sterile neutrinos and the light neutrinos, and the 
elements |J7e4| and |J7e5| of the leptonic mixing matrix, 
which describe the mixing of the electron neutrino flavor 
with the heavy neutrino mass states 1/4 and 1/5. Obvi- 
ously, for the 3-1-1 case, only 1/4 is present. The best fit 
points for the two scenarios are summarized in Table [H 
To illustrate the impact of sterile neutrinos on the fit 
to reactor anti-neutrino data graphically, we compare in 
Fig. [T] the data to the predictions for the no oscillation 
case (green) and the best fit 3+1 (blue) and 3-1-2 (red) 
models. Note that, even for no oscillations, the prediction 
may deviate from 1 due to nuisance parameters included 
in the fit, parameterizing systematic uncertainties. The 
fit is dominated by Bugey3 spectral data at 15 m and 
40 m and the precise rate measurement from Bugey4. 

In the lower part of Fig. [2] we show the of the SBL 
reactor fit as a function of Am^j^ . Using the new flux pre- 
dictions (solid curves) we flnd a clear preference for sterile 
neutrino oscillations: the Ax^ between the no oscillation 
hypothesis and the 3-1-1 best fit point is 8.5, which im- 
plies that the no oscillation case is disfavored at about 
98.6% CL (2 dof). In the 3+2 case the no oscillation hy- 
pothesis is disfavored compared to the 3+2 best fit point 
with Ax^ = 12.1, or 98.3% CL (4 dof). In contrast, with 
previous fiux predictions (dashed curves) the improve- 
ment of the fit is not significant, with a Ax^ between the 
best fit points and the no oscillation case of only 3.3 and 
5.4 for the 3+1 and 3+2 hypotheses, respectively. 

Global analysis of SBL data. The constraints from 
the reactor experiments under discussion play an impor- 
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Figure 2: x^ from global SBL data (upper panel) and from 
SBL reactor data alone (lower panel) for the 3+1 (blue) and 
3+2 (red) scenarios. Dashed curves were computed using the 
old reactor anti-neutrino flux prediction Q, sohd curves are 
for the new one [s] . All undisplayed parameters are minimized 
over. The total number of data points is 137 (84) for the global 
(reactor) analysis. 



tant role in a combined analysis of all SBL oscillation 
data, including the LSND and MiniBooNE anomalies. 



LSND has provided evidence for t'^ 



i>e transitions 0], 
and MiniBooNE has reported an excess of events in the 
same channel, consistent with the LSND signal This 
hint for oscillations is however not confirmed by a Mini- 
BooNE search in the i^^ — > i^e channel [l^, where the 
data in the energy range sensitive to oscillations is con- 
sistent with the background expectation. These results 
seem to suggest an explanation involving CP violation 
in order to reconcile different results for neutrino and 
anti-neutrino searches. An explanation of the LSND and 
MiniBooNE anomalies via sterile neutrino oscillations re- 
quires the mixing matrix elements |J7e4| and/or |J7e5| to 
be non-zero. Reactor experiments are sensitive to these 
parameters, and while analyses using previous flux pre- 
dictions lead to tight constraints on them, the new fluxes 
imply non-zero best fit values (Tab.|l|) and closed allowed 
regions at 98% CL. Hence, the interesting question arises 
whether a consistent description of the global data on 
SBL oscillations (including LSND /MiniBooNE) becomes 
now possible. To answer this question we perform a fit 
by including, in addition to the reactor searches for De 
disappearance, the LSND ^ and MiniBooNE 0, [H] re- 
sults, as well as additional constraints from ap pea rance 
experiments 2^, 21 1, j/^ disappearance searches |22|. and 
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Figure 3: Global constraints on sterile neutrinos in the 3+1 
model. We show the allowed regions at 90% and 99% CL from 
a combined analysis of the LSND 0] and MiniBooNE anti- 
neutrino [3 signals (filled regions), as well as the constraints 
from the null resuhs of KARMEN NOMAD [H and 

MiniBooNE neutrino 19] appearance searches (blue contour). 
The limit from disappearance experiments (green contours) 
includes data from CDHS [22], atmospheric neutrinos, and 
from the SBL reactor experiments. For the latter we compare 
the results for the new anti-neutrino flux prediction from 
(solid) and the previous ones (dashed). The region to the 
right of the curves is excluded at 99% CL. 



atmospheric neutrinos. Technical details of our analysis 
can be found in [1, [l3l and references therein. 

In the 3-fl scheme the SBL experiments depend on 
the three parameters Amli, \Ue4\, and \U^4\. Since 
only one mass-scale is relevant in this case it is not 
possible to obtain CP violation. Therefore, oscillations 
involving one sterile neutrino are not capable of rec- 
onciling the different results for neutrino (MiniBooNE) 
and anti-neutrino (LSND and MiniBooNE) appearance 
searches. Fig.[3]compares the allowed regions from LSND 
and MiniBooNE anti-neutrino data to the constraints 
from the other experiments in the 3-1-1 model. Note 
that, even though reactor analyses using the new flux 
prediction prefer non-zero Ud, no closed regions ap- 
pear for the disappearance bound (solid curve), since 
sin^20SBL = 4]J7e4p|C^At4p can still become zero if 
Uf_i4 = 0. We find that the parameter region favored by 
LSND and MiniBooNE anti-neutrino data is ruled out by 
other experiments, except for a tiny overlap of the three 

1 eV^. Note that in 
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99% CL contours around Af 
this region the constraint from disappearance data does 
not change significantly due to the new reactor flux pre- 
dictions. Using the PG test from [2^ we find a compat- 
ibility of the LSND-KMiniBooNE(i^) signal with the rest 
of the data only of about 10-^ with Xpg = 21.5(24.2) 



S+2 0.47 0.128 0.165 0.87 0.138 0.148 1.64 110.1/130 
1+3+1 0.47 0.129 0.154 0.87 0.142 0.163 0.35 106.1/130 

Table II: Parameter values and the global best fit 

points for 3+2 and 1+3+1 oscillations (Am^'s in eV^). 
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Figure 4: Predicted spectra for MiniBooNE data and the 
transition probability for LSND (inset). Solid histograms re- 
fer to the 3+2 global best fit point (Tab. |II]), dashed his- 
tograms correspond to the best fit of appearance data only 
(LSND, MiniBooNE u/u, KARMEN, NOMAD). For Mini- 
BooNE we fit only data above 475 MeV. 



for new (old) reactor fluxes. Hence we conclude that the 
3+1 scenario does not provide a satisfactory description 
of the data despite the new hint coming from reactors. 

Let us move now to the 3+2 model, where SBL exper- 
iments depend on the seven parameters listed in Tab. [Ill 
In addition to the two mass-squared differences and the 
moduli of the mixing matrix elements, also a physical 
complex phase enters, S = a,rg{Uf_i4U*^U*^Ue5)- This 
phase leads to CP violation in SBL oscillations [1, , 
allowing to reconcile differing neutrino and anti-neutrino 
results from MiniBooNE/LSND. Tab. [Ushows the para- 
meter values at the global best fit point and the corre- 
sponding value. Changing from the previous to the 
new reactor flux calculations the decreases by 10.6 
units, indicating a significant improvement of the descrip- 
tion of the data, see also upper panel of Fig. [2] From that 
figure follows also that going from 3+1 to 3+2 leads to 
a significant improvement of the fit with the new reactor 
fluxes, which was not the case with the old ones. The 
improves by 11.2 units, which means that 3+1 is dis- 
favoured at the 97.6% CL (4 dof) with respect to 3+2, 
compared to Ax^ = 6.3 (82% CL) for old fluxes. 

In Fig. [1] we show the prediction for the Bugey spectra 
at the global best fit point as dashed curves. Clearly they 
are very similar to the best fit of reactor data only. Fig. [4| 
shows the predicted spectra for MiniBooNE neutrino and 
anti-neutrino data, as well as the LSND -+ transi- 
tion probability. Again we find an acceptable fit to the 
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LSND+MB(!/) vs rest appearance vs disapp. 



old 


new 


old 


new 


xl>G,3+2/dof 25.1/5 


19.9/5 


19.9/4 


14.7/4 


PG3+2 10"* 


0.13% 


5 X lO"'* 


0.53% 


xi.G,i+3+i/dof 19.6/5 


16.0/5 


14.4/4 


10.6/4 


PGi+3+1 0.14% 


0.7% 


0.6% 


3% 



Table III: Compatibility of data sets [23| for 3+2 and 1+3+1 
oscillations using old and new reactor fluxes. 



data, although in this case the fit is slightly worse than 
a fit to appearance data only (dashed histograms) . Note 
that MiniBooNE observes an event excess in the lower 
part of the spectrum. This excess can be explained if only 
appearance data are considered, but not in the global 
analysis including disappearance searches 0. Therefore, 
we follow [Toj and assume an alternative explanation for 
this excess, e.g. In Tab. |TTT]we show the compat- 

ibility of the LSND/MiniBooNE(j>) signal with the rest 
of the data, as well as the compatibility of appearance 
and disappearance searches using the PG test from (23j . 
Although the compatibility improves drastically when 
changing from old to new reactor fluxes, the PG is still 
below 1% for 3+2. This indicates that some tension be- 
tween data sets remains. We considered also a "1+3+1" 
scenario, in which one of the sterile mass eigenstates is 
lighter than the three active ones and the other is heav- 
ier [2^. As can be seen from Tabs. and lllll the fit 
of 1+3+1 is slightly better than 3+2, with Ax^ = 15.2 
between 3+1 and 1+3+1 (99.6% CL for 4 dof). How- 
ever, due to the larger total mass in neutrinos, a 1+3+1 
ordering might be in more tension with cosmology than 
a 3+2 scheme (27l - [29| . Fig. [5] shows the allowed regions 
for the two eV-scale mass-squared differences for the 3+2 
and 1+3+1 schemes. 

Discussion. Let us comment briefly on other signatures 
of eV sterile neutrinos. We have checked the fit of solar 
neutrino data and the KamLAND reactor experiment, 
and found excellent agreement. The effect of non-zero 
values of 11^4 and Ue5 for these data are similar to the 
one of J7e3 in the standard three-active neutrino case, and 
hence the 3+2 best fit point mimics a non-zero LL3 close 
to the preferred value of these data, see [l|, 0, • The 
MINOS long-baseline experiment has performed a search 
for sterile neutrinos via neutral current (NC) measure- 
ments [31}. We have estimated that the best fit points 
reported in Tab. [Til lead to an increase of the of MINOS 
NC data as well as Xpg ^'Y a few units [s^ . Radioactive 
source measurements in Gallium solar neutrino experi- 
ments report an event deficit which could be a manifes- 
tation of electron neutrino disappearance due to eV-scale 
sterile neutrinos [sl] . Our best fit points fall in the range 
of parameter values found in |32j capable to explain these 
data. Finally, eV-scalc sterile neutrinos may manifest 
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Figure 5: The globally preferred regions for the neutrino 
mass squared differences Am^i and Am§i in the 3+2 (upper 
left) and 1+3+1 (lower right) scenarios. 



themselves in cosmology. Recent studies [27h29| indicate 
a slight preference for extra radiation content in the uni- 
verse (mainly from CMB measurements) and one or two 
sterile neutrino species with masses in the sub-eV range 
might be acceptable. Big-Bang nucleosynthesis leads to 
an upper bound on the number of extra neutrino species 
of 1.2 at 95% CL [s^l, which may be a challenge for two- 
sterile neutrino schemes, or indicate a deviation from the 
standard cosmological picture. 

In conclusion, we have shown that a global fit to short- 
baseline oscillation searches assuming two sterile neutri- 
nos improves significantly when new predictions for the 
reactor neutrino flux are taken into account, although 
some tension remains in the fit. We are thus facing an 
intriguing accumulation of hints for the existence of ster- 
ile neutrinos at the eV scale, and a confirmation of these 
hints in the future would certainly be considered a major 
discovery. 
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